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We consider a model of the state evolution of relativistic vector bosons, which includes both the 
dynamical equations for the particle four-velocity and the equations for the polarization four-vector 
evolution in the field of a nonlinear plane gravitational wave. In addition to the gravitational 
minimal coupling, tidal forces linear in curvature tensor are suggested to drive the particle state 
evolution. The exact solutions of the evolutionary equations are obtained. Birefringence and tidal 



Q ' deviations from the geodesic motion are discussed. 

in ; 

: 1 INTRODUCTION 

(N . 
> ■ 

^ I The formulation of the covariant dynamic equations for the particles with internal structure 

O ■ OT with supplementary degrees of freedom unavoidably involves consideration the tidal forces, 

O ■ 

^ ■ i. e., forces linear in the curvature tensor. 

^ I A. Papapetrou in the paper'^^\ deriving the covariant equations for the spinning particle, 

O . presented the first example of the tidal forces. Considering the multipole representation 

^ • of the interaction between the particle possessing internal structure and the external field, 

W. G. Dixon*^^) and A. H. Taub*^^-* introduced the coupling of the curvature and the particle 
■ quadrupole moment. W. Israel'-'^-* formulated the covariant dynamics of the macroscopic 

;h ' polarization in the medium also using the tidal interactions. At present the covariant theory 

of the dynamics of spinning objects is well developed due to the interest to the problem of 
emission of the gravitational waves (e.g., Ref. 5). 

As it was shown by I. T. Drummond and S. J. Hathrell in Ref. 6, the tidal terms ap- 
pear in the covariant electrodynamics. The birefringence induced by curvature^^"*) is an 
electrodynamical phenomenon, which also admits interpretation in terms of the photon mo- 
tion under the influence of tidal forces. This interpretation is based on the fact that due 
to gravitationally induced birefringence the energy and the instant direction of the photon 
momentum vector depend on the direction of the polarization vector, rotating in the gravity 
field. 

On the other hand, when the electromagnetic fields dominate and the tidal forces are 
negligible in comparison with them, we can use the relativistic Bargmann-Michel-Telegdi 
equations*^^-* , describing the minimal interaction with gravitation and do not displaying the 
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tidal forces. Nevertheless, as it have been emphasized by I. B. Khriplovich*^^°\ there ex- 
ists a direct analogy between electromagnetic and tidal coupling of the particle spin. One 
can introduce, as was advocated by Khriplovich*^^°\ new part of the Hamiltonian, namely, 
— jRabcdS"''^S^'^, where Rated is a Riemann curvature tensor and is a spin tensor, instead of 
(or in addition to) the electromagnetic coupling part ^F^bS'^^. In other words, both electro- 
magnetic and tidal coupling of spin can be described analogously. Finally, I. B. Khriplovich 
noted in Ref. 10 that the developed formalism and the corresponding equations are applicable 
to the vector bosons with the spin equal to one. 

We consider this idea as an initial point for our investigations and we use the evolutionary 
equations for vector boson with vanishing electromagnetic field and dominating tidal forces 
(see the Section 2 of this paper) . In the Section 3 we integrate the dynamic equations exactly 
for the case of nonlinear gravitational wave background and present two particular simplified 
exact toy-models. 

2 EVOLUTIONARY EQUATIONS 

The covariant model of the vector boson evolution contains the pair of equations. 

DJji 1 

= -—R%iJJ'^'u-^. (1) 

i?*?,,^H^S'C/-. (2) 
Ut mc 

The first equation is a dynamical one for the determination of the particle four- velocity time- 
hke vector = ^ {U'^Ui = 1). The second equation is the equation of evolution of the 
polarization space-like four- vector (H*Sj = — -E^, the constant E has the dimensionality 
of the Planck constant). D denotes the covariant differential, r is the affine parameter along 
the particle world-line (particle proper time). The tensor R*%i^ is a right-dual to Riemann 
tensor 



Dr mc 

DE' 1 



-klm — 2 -fep? ••'wi' 



The Levi-Civita tensor e^^f^ is defined as usual with 

— ^ ^ E^'^ , ^pqlm — y/ 9 ^pqlmt (4) 

where E^'^'-"^ is a completely skew-symmetric symbol with E^^"^^ = — -E0123 = 1- 

The equation (1) looks like the relativistic dynamic equation mc^^ = with the force 
four- vector expressed in terms of the Riemann tensor. It is a typical "tidal force". This 
force is evidently orthogonal to the four- velocity vector (F'[/j = 0), and this fact combined 
with the velocity normalization law admits us to speak about particle mass conservation. 



3 



Analogously, one can consider the equation (2) as an equation with "tidal force" acting 
on the polarization four-vector. The right part of this equation is apparently orthogonal to 
the S* vector. Consequently, S'^Sfc = —E"^ is a constant of motion. 

Finally, one can see that due to (1) and (2) 

Dr^ ^ Dr Dr ^ ^ 

i.e. the scalar product (f/*Sj) = const is an integral of motion. This constant is considered 
to be equal to zero in order to display the law of orthogonality of the four- velocity vector and 
polarization four- vector for the massive vector boson as well as for the massless photon*^^°^. 

One can note that S** and quantities are considered to be pseudo- vectors, nevertheless, 
the product of the right-dual Riemann tensor with spin or polarization is a true tensor. For 
the sake of simplicity we shall omit the prefix "pseudo" below. 

2.1 Remark on the derivation of the evolutionary equations 

The simplest way to obtain the evolutionary equations (1), (2) is the following. Let us 
consider the Bargmann-Michel-Telegdi equations*^^) with gyromagnetic ratio g — 2, i.e. in 
the case of the absence of the anomalus magnetic moment: 

^ = ^nu\ ^ = ^,FiS\ (6) 
Dt mc^ Dt mc^ 

Then, by the analogy, discussed by I. B. Khriplovich^^"^ , let us replace the term describing 
the contribution of the electromagnetic field :^FJ'^, by the tidal contribution —-^^R^.j^i^S^^ . 
If we express the spin-tensor 5"'"* in terms of spin vector 

gin. ^ ^Inmu^S,, = ^Qfc^^C/'-S'"^, (7) 

use the definition of the right-dual Riemann tensor (3), and then replace the spin four- vector 
by the polarization four-vector S*, we shall obtain the equations (1) and (2). 

2.2 Gravitational wave background 

The model of the plane gravitational wave (GW) background based on the exact solution 
of the Einstein equations in vacuum^^^) is well-known and highly fruitful. The metric of the 
plane GW can be represented in the following form: 

ds" = 2dudv - L{uf{[e^^^''\dx'f + e-^'^^''\dx''f] cosh2-f{u) + 2 sinh2-t {u)dx^dx''), (8) 

where 



u 



= {ct - x^)/V2, v={ct + x^) /V2 (9) 
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are the retarded and the advanced time, respectively. The functions L(ti), 7(ti) are 
coupled by one equation only*^^^) 

L"{u) + L{u) {[(3\u)fcos}i^^{u) + W{u)f} = 0, (10) 

which is the unique remaining nontrivial Einstein's equation for the case of vacuum. The 
prime denotes the derivative by the retarded time u. The functions (5{u) and ^{u) are 
assumed to be arbitrary ones. We shall say that the gravitational wave is of the first polar- 
ization, if I3{u) 7^ 0, ^{u) — 0, and of the second polarization, if I3{u) — 0, ^{u) ^ 0. 
The initial conditions for the functions L{u), I3{u),^{u) are the following: 

L(0) = 1, L'(0) = 0, /?(0)=7(0)=0, /3'(0) = 7'(0) = 0. (11) 

The Riemann tensor calculated with the metric (8) has only four nontrivial components 

-Rl^^ = i?^„3„ = L-2[L2/3'cosh2 27]', 

Rl^u = -L-^iL'e^^ii - /3'sinh27COsh27)]', 

Rl^u = -L-^[L2e-^^(7' + /3'sinh27COsh27)]'. (12) 
The components of the right-dual Riemann tensor are equal, correspondingly, to 

r>* 7" 2 p3 p* 7" 2 p2 

^^2u2u ~ ^ ' -^■u2u^ -^3u3u ~ ^ ' -^-uSu' 

p* p* 7" 2 p2 t2 p3 /'I 

-"-2m3m — ^3u2u — ^ ' ^■u2u — ^ ' ^■u3u ■ \^^) 

The space-time with metric (8) admits the existence of five Killing's vectors, three of them, 
namely, 

C{v) = C{2) = ^2: C(3) = (14) 

form the Abelian subgroup. The vector ^^^^ possesses the following properties: 

i.e., it is a null covariantly constant vector, orthogonal to the Riemann tensor and to its dual 
tensors. 

3 THE EXACT SOLUTION FOR THE MODEL DYNAMICAL SYSTEM 
3.1 The first integrals of motion 

We have found (see Section 2) that due to the structure of the evolutionary equations three 
quadratic quantities 

U'Ui = const = 1, E0 = const = -E'^, EiW = const = (16) 
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are constant along the particle world-line independently on the gravitational background 
properties. 

Next two integrals of motion 

Uv Qik C{v)U^ = const = Cy, Ey == C{v)-^ = const = Ey (17) 

exist due to the GW symmetry. The convolution of the equations (1) and (2) with QikC'^y-^ 
gives the formulae (17), if we use the properties (15) of the null covariantly constant Killing 
vector. 

Finally, accounting that 

^ = = = a, (18) 
ar 

we can link the particle proper time r with the retarded time u: 

u 

T = — + const. (19) 
Cy 

Note, that for massive boson 0. We shall use the relationship (19) with the con- 

stant equal to zero in order to reparametrize the remaining differential equations for the 
components of [/' and S', parallel to the GW front plane. 



3.2 The motion in the GW front plane 

Let us extract four equations from (1) and (2) by means of convolution of (1),(2) with 
the second and the third KiUing vectors from the abehan subgroup (14). Denoting the 
corresponding convolutions by 

Ua = gik^Uu\ Ea = 9ikCi.)'^\ (20) 

where the Greek indices run from 2 to 3, we obtain the following system of coupled equations: 

^ Ku,u9AC.Mu) - EyUp{u)l (21) 



1 E 

I ( \ TD* 

2 5*07 1""^ ~ ^au-iu 



g'^^[CyE0{u)-EyUp{u)]. (22) 



The evident symmetry of the equations (21), (22) allows us to introduce some new unknown 
functions 

X^{u) = CyE^iu) - EyU^{u). (23) 

Then we obtain the following two-dimensional key subsystem as the differential consequence 
of (21),(22): 

' X'A^(al «| \ ( X2 
X',] \alal]\x^ 



(24) 
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Here 

al {u)^ — + p' cosh^ 27, al {u) = e^^ [7' - sinh 2j cosh 2j] , 
Ij 

u 

al{u)=e-'^f^ [7' + /3' sinh 27 cosh 27] , a\{u) = — - cosh^ 27. (25) 

It is interesting to mention that the key subsystem happens to be self-closed and does not 
contain the tidal terms with the Riemann tensor. Such two-dimensional subsystems were 
considered, solved and used in Refs. 8 and 12. The solution of (24), (25) is the following 

X2 \ ( I I sinh 7 I / cos^ — sin^ \ / X2(0) \ ^^^^ 

J y J y sinh 7 cosh 7 J y sin^^ cosV' J y -^3(0) J 

where the function ip{u) has the form 

u 

iP{u) = J (3' smh2-fdu. (27) 


The initial data Xa{0) are predetermined by those for the four- velocity and polarization 
vectors: 

Xa{0) = EaCy — E^Ca, Ca = UaiO), Ea = 'B.a{0). (28) 

Then the solution of the total four-dimensional subsystem of equations (21), (22) can be 
reduced to the quadratures 

^a{u) = ^ [XM + E,UM] , (29) 



Ua{u) ^C^-— [du Rl^^M9'^{n)X0{u). (30) 



In (8), (12), (13) and (26) one can find all the functions which are necessary for the integration 
in (30). 

3.3 The quadratures for the remaining unknown functions 

Now, when X„ are obtained and Ua{u) as well as ^^(m) are represented via X^, we can 
extract Uu and X^ from the quadratic integrals (16): 

Uu{u) = - g''''{u)U^{u)U,{u)l a 7^ 0, (31) 

^u{u) = + g'"'{u)'EM^Ml ^ 0. (32) 

Finally, the particle position as a function of the retarded time is described by the formulae 

u 

x"(k) = x"(0) + ^ 1 g''''{u)U0{u)du, (33) 
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1 " 

v{u) = -^(0) + ^ / Uu{u)du. 



(34) 



Thus, the formulae (17),(19),(29)-(34) represent in the quadratures the exact solution of the 
equations (1),(2), containing 12 functions of the proper time r. 



3.4 The first exact toy- model 

If the massive boson was at rest at the initial moment ti = in the chosen reference frame, 
i.e. 

^2 = = 0, = ^ = t^«(0), (35) 

and if the GW has only the first polarization, i.e. ^{u) — 0, then ■?/;(«) = and we can 
simplify the formulae (26), (30), (29): 



S,(«) = Lel^E, - ^ (Le«)' , E,{u) = Le'^E, + ^ (Le-«)' . 

UlC ^ ' 'mr ^ ' 

The relationships (16) yield 



El 



mc 



(36) 

(37) 
(38) 



2E^ + Ei + Ei = E^, (39) 

i.e. one of the constants E^, E2 or E3 can be expressed in terms of other ones and of the E 
normalization constant. 

Computing the particle energy 



mc 



£(u) = mc'Uo = ^(C/„ + U^) 



mc 



V2 



2\/2Cy 



[l-g^^UJJp + 2Cl), 



(40) 



and using (37) we obtain 



Eiu) — mc? — - — 
4m 



(41) 



Finally, let us calculate the Ui component of the four-velocity vector: 

2 



(42) 
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3.5 The second toy-model 



Let the GW be of the first polarization as well. Let at the moment u — the polarization 
three-vector of a particle be directed along the Ox^ axis, and the initial particle motion be 
one-dimensional, for example along Ox^ axis. It is possible, if the constants are chosen as 
follows: 



Then we see that 



U2{u) = C2, Us{u) 



C2E 



Ui{u) = 



EC2 



1 - 



\/2mc 

2/3 \ ^2 /^/ 



L2 



2m? \ L 



EU. 



The formulae of two last subsections will simplify the next discussion. 



(43) 

(44) 
(45) 

(46) 



DISCUSSION 



The obtained exact solution of the evolutionary model demonstrates explicitly the following 
general properties. 

(i) Birefringence induced by curvature 

Using the formulae (40) with (30) or an explicit formula (42) we can conclude that the 
particle energy at the moment u depends on the initial data for the polarization four- vector. 
Since the polarization four-vector is normalized and orthogonal to the four-velocity (see 
(16)), only two initial components of this four- vector are independent (see (39)). In general 
case it is convenient to choose the E2 and £'3 parameters to be independent ones. Since 
the particle energy depends on two initial polarization parameters, we deal with the bosonic 
analogue of the optical birefringence. Since the polarization parameters are involved into 
the particle energy formula just due to the tidal interaction, we can denote the effect of this 
type as birefringence induced by curvature^^"^'' . We can compute the energy shift (41) for 
i?2 = and for E^ = Q. The corresponding ratio of the polarizationally induced energy shifts 
is equal to 

[£ - mc%,=o : [£ - mc\E,=o = {L' + L(5'f : (L' - L(5'f. (47) 

Note that the polarizationally induced energy shifts for the first toy-model are of the same 
sign. In the weak GW field they differ only in the second order of (3. Nevertheless, in the 
case of initially moving particle these shifts differ more considerably and have contributions 
linear in /3. 
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The second toy-model demonstrates the degenerate case: if the E2 and components 
vanish, i.e. the initial direction of the polarization three-vector coincides with that of the 
GW propagation, then the energy depends only on the polarization four-vector modulus, 
and birefringence is not displayed. 

(ii) Particle nongeodesic motion 

In the presence of the tidal interaction the particle motion is no longer the geodesic one. 
The tidal force in the right part of the formula (1) produces the particle acceleration and 
the rotation of the four- velocity vector. The formulae (30), (31), (26) describe explicitly these 
phenomena. Speaking about the particle rotation, we emphasize the following feature. If 
the initial direction of the particle motion is along the Ox^ axis (see the second toy - model) , 
then in the field of tidal forces the U2 component remains constant, but the U^lu) and Ui{u) 
components appear. Certainly, it is not a standard rotation, because the instantaneous 
radius of the quasi-orbit grows. But one can find some analogy with the rotation of the 
charged particle in the nonhomogeneous magnetic field. 

In principle, the phenomena of the particle nongeodesic acceleration and rotation are 
the supplementary contribution to the geodesic ones, caused by the minimal coupling with 
gravitation. This statement can be illustrated by the formula (45), in which the term in the 
parentheses describes the geodesic variation of the longitudinal velocity, and the last term 
appears due to the tidal interactions. 

(iii) Rotation of the polarization vector 

The effect of the rotation of the polarization vector is displayed by the formulae (29) 
and (26). The third two-dimensional matrix in the product in the right part of (26) is the 
standard rotation matrix, the ip{u) being the phase of rotation. Even if we shall neglect 
the term with right-dual Riemann curvature tensor in (30), nevertheless, the polarization 
four- vector (see (29)) will rotate due to the properties of the Xa{u) (26). In addition, the 
rotation of the particle three- velocity produces the supplementary rotation of the polarization 
three-vector, which is induced originally by curvature. 

Thus, we can conclude that the gravitational waves effect on the vector boson can be 
considered as a combination of geodesic and tidal acceleration and rotation, leading to the 
birefringence phenomenon. 
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